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O \ ABSTRACT 

We have detected the CO(2-l) transition from the submillimetre galaxy (SMG) 

^ ■ LESS J033229.4-275619 at z = 4.755 using the new Compact Array Broadband Backend 

\ system on the Australian Telescope Compact Array. These data have identified a massive 

• *h . gas reservoir available for star formation for the first time in an SMG at z ~ 5. We use the 

' luminosity and velocity width (FWHM of ~ 160 km s^ 1 ) of the CO(2-l) line emission to 

constrain the gas and dynamical mass of M gas ~ 1.6 x 1O 1O M0 and Md yn (< 2kpc) ~ 
5 x 10 10 (O.25/sin 2 i)M , respectively, similar to that observed for SMGs at lower redshifts 
of z <~ 2^1, although we note that our observed CO FWHM is a factor of ~ 3 narrower than 
typically seen in SMGs. Together with the stellar mass we estimate a total baryonic mass of 
Mbary — 1 x 10 11 M , consistent with the dynamical mass for this young galaxy within the 
uncertainties. Dynamical and baryonic mass limits of high-redshift galaxies are useful tests of 
galaxy formation models: using the known z ~ 4-5 SMGs as examples of massive baryonic 
systems, we find that their space density is consistent with that predicted by current galaxy 
formation models. In addition, these observations have helped to confirm that z ~ 4-5 SMGs 
possess the baryonic masses and gas consumption timescales necessary to be the progenitors 
of the luminous old red galaxies seen at z ~ 3. Our results provide a preview of the science 
that ALMA will enable on the formation and evolution of the earliest massive galaxies in the 
Universe. 

Key words: galaxies: high-redshift - galaxies: evolution - galaxies: formation - submillime- 
tre - galaxies: individual LESS J033229.4-275619 



1 INTRODUCTION 

Blank field surveys at millimeter and submillimetre wavelengths 
have uncovered a population of extremely luminous, highly op- 
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tically obscured star-f orming galaxies, so termed submillimetre 
galax i es (SMGs; e.g. ISmail. Ivison & BlairJ 1997 ; iHughes et al.l 
1 19981 : iBarger et al.1 Il998l : ICo ppin et alj 120061 : IWeifi et al.l l2009h . 



SMGs are infrared luminous (Lir ^ 5 x 10 L©) and have high in- 
ferred star formation rates (SFRs) in excess of lOOOM yr"\ and 
the redshift distribution of radio-identified 850-/im selected SMGs 
(repre senting ~ 60% of t he SMG population; e.g. Ilvison et"al] 
120071) peaks at z ~ 2.2 l lChapman etal] 120051). There are in- 
dications that SMGs are strongly clustered dWebb et alj 120031 ; 
iBlain et alj |2004 IWeiB et al.l l2009h . while molecular gas obser- 
vations of SMGs have demonstrated that they contain the mas- 
sive compact gas reservoirs needed to form a lum inous spheroid 
dGreve et al.ll2005l : lTacconi et alj2006ll2008ll20irj|) . The combina- 
tion of these properties make them attractive candidate s to be pro- 
genitors of today's most massive elliptical galax ies (e.g. lLillv et al.l 
Il999t iGenzel et alj|2003l ; ISwinbank et aT]|2006h . 

While the bulk of SMGs lie at z ~ 2, recently the first exam- 
ples of submm- selected SMGs at z ~ 4-5 have been uncover ed 
ilCapaketalJ [20081: ISchinnerer etal] 120081 ; iDaddi et alj l2009aMbl 
Copp in et alj 120091 : llCnudsen et al.1 |2010|) . reviving the potential 



worry that significant numbers of SMGs could li e at z > 4 and 
are missed from curren t spectroscopic surveys (e.g. lChapman et all 
120051 : IPope e7ai]|2005h . The presence of even a small number of 
galaxies at z > 4 forming stars at a rate of ~ lO 3 M0yr _1 
could pose a serious problem for current hierarchical models, 
which predict a z > 4 SMG surface density which is barely 
consistent with highl y incomplete cu rrent observational estimates 
(~ 10deg~ 2 ; see lCoppin et al.l 12009). However, the galaxy lumi- 
nosities are not robust predictions of the models since they depend 
upon t he adopted initial m ass function (IMF) and dust properties 
(e.g. lSwinbank et al1 l2008). Much stronger and more reliable con- 
straints on the models are provided by comparisons of the bary- 
onic (gas+stars) and dynamical masses fo r these young galaxies 
(e.g. lGenzel et alfeOOalBaugh et all2"005l) . 

In this letter we report on the detection of CO(2-l) emission 
from a luminous SMG, LESS J033229.4-275619 (LESS J033229.4 
hereafter), discovered in the deep 870-/im Large Apex Bolome- 
ter CAmera (LABOCA) Extended Chand ra Deep Field South 
(ECDFS) Survey (LESS; IWeifi et alj|2009t) . LESS J033229.4 has 
a deboosted 870-/im flux density of 5.1 ± 1.4mJy (4.6a sig- 
nificance), a faint radio and mid-infrared counterpart, and opti- 
cal spectroscopic and photometric observations identif ying it as an 
ultraluminous far-infrared galaxy at z — 4.76 (see ICoppin et al.1 
2009). This is currently the highest redshift SMG known, suggest- 
ing that intense star formation was occurrin g in galaxies when th e 
Universe was just over 1 Gyr old (see also iRiechers et al]|2009l) . 
Here, we use observations of CO(2-l) to a) confirm the counter- 
part identification and redshift of the source; b) to investigate the 
Lco-Lir relation for SMGs at z ~ 5; c) to determine if z > 4 
SMGs have the baryonic and gas consumption time scales neces- 
sary to be the progenitors of the lu minous red galaxies found at 
z ~ 3 jMarchesini etai]|200l l2009); and d) to derive a total bary- 
onic mass for the system by c ombining ou r stel lar and gas mass 
estimates, in order to test the iBaugh et al.1 d2005t) galaxy forma- 
tion model. We adopt cosmological parameters of SIa = 0.73, 
fi m = 0.27, and H = 71kms _1 Mpc" 1 . At z = 4.76 the an- 
gular scale is 6.5 kpc arcsec -1 , and the Universe is 1.3 Gyr old. 

2 OBSERVATIONS AND DATA REDUCTION 



—27° 56 19.40 ) was observed in 2009 September over a period 
of 6 nights with the Australia Telescope Compact Array (ATCA), 
fitted out with the recently upgraded Com pact Array Broadband 
Backend (CABB; F ems & Wilson 2002). For improved (u, v)- 
coverage, we employed a combination of the two most sensitive 
and compact 5-antenna hybrid configurations, H75 and H168. We 
utilised the maximum available CABB 2.048 GHz bandwidth (en- 
suring that the CO line of interest would lie well within the spectral 
window), with 2048 channels across the intermediate-frequency 
(IF) band, corresponding to a 1 MHz spectral resolutio n. Assum- 
ing a source redshift of z = 4.762 jCoppin et al1 l2009). we tuned 
the 7 mm IF1 receiver to a central frequency of 40.010 GHz, re- 
sulting in a frequency coverage of 39.035-^10.984 GHz, covering 
CO(2-l) emission between z = 4.625-4.906. The system tem- 
perature ranged from 90-160 K over the observing run. Phase, am- 
plitude, and bandpass calibration were acquired by short observa- 
tions of PKS 0402-362 every 10 mins, with pointing performed 
on the same source every hour. Flux density calibration was per- 
formed using short observations of Uranus each night (and Mars in 
one instance), the recommended standard ATCA calibrator source, 
and we estimate the calibration uncertainty to be <; 20% based 
on the currently avail able standard MIRIAD model for Uranus 
dSault&Killeenll999l : see lEmonts et all20ld) . 

The data w ere calibrated, mapped and analysed us ing the stan- 
dard MlRIArfl JSault & Kill een 1999) and KARM/J 2 ! jGooch et al.1 
1 19961) packages. The corresponding synthesized beam adopting a 
natural weighting is ~ 12x9 arcsec . Overall 35 hr of observations 
were obtained, including 25 hr of on-source integration time. The 
resultant 1 a rms channel noise in the 1 MHz-wide channel spec- 
trum (corresponding to ~ 7kms~ 1 ) is ~ 0.44mJybeam _1 , which 
is consistent with sensitivities reached by other 7 mm ATCA/CABB 
experiments (e.g. lEmonts et al.ll2010j) . 

3 OBSERVED MOLECULAR GAS PROPERTIES 

Inspection of the combined data cube reveals the presence of a 
relatively narrow emission line near the phase centre, blueshifted 
by ~ 400kms _1 from the central frequency. To test the reality 
of the detection, we split the visibility data into two halves and 
find that the line persists in the data cubes at lower signal-to-noise. 
To quantify the significance and parameters of the CO(2-l) line 
emission, the visibilities were resampled to a velocity resolution of 
90 km s -1 (chosen to yield optimal line SNR), providing a 1 a rms 
channel sensitivity of ~ O.MmJybeam" 1 , revealing a ~ 5a de- 
tection (see Fig.[TJ, the brightest peak in the entire data cube with 
this binning. In Fig. [2] we plot contours of the emission for the line 
at the brightest 90kms— 1 channel overlaid on a true-colour near- 
infrared HAWK-I JK image of LESS J033229.4 from Zibetti et al. 
(in preparation). A Gaussian profile fit to the CO spectrum indicates 
a full width at half maximum (FWHM) of 160 ± 65 km s" 1 , and a 
velocity integrated line flux of Fco(2-i) = 0.09 ± 0.02 Jy km s~ . 
The uncertainties on the CO line parameters have been estimated 
by rebinning the data cube to several different spectral channel res- 
olutions from 50-200 km s _1 and finding the 1 a spread of each 
parameter, which are comparable to the individual parameter fit 
uncertainties. We determine a best-fitting CO(2-l) redshift of z = 
4.755±0.001, corresponding to an offset of -370±30km s" 1 w ith 
respect to the optical redshift derived by IVanzella et alj d2006t) of 
z = 4.762±0.002 determined from the asymmetric Lya emission, 
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1 http://www.atnf.csiro.au/computing/software/miriad 
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Figure 1. CO(2-l) spectrum of LESS J033229.4 binned into 90kms _1 
channels, extracted through the brightest pixel from Fig. [5] showing a 
~ 5 a line detection overlaid with a best-fitting Gaussian with a FWHM 
of 160 ±65 kms -1 . The upper axis indicates the redshift of the CO(2— 1) 
transition, from which the line emission fit indicates a z = 4.755 ± 0.001. 
The CO emission line is blueshifted by fa 370 kms -1 from the optical 
Lya redshift (dotted vertical line and error bar) and consistent with a red- 
shift of z = 4.751 ± 0.005 derived from a VLT [Oil] emission spectrum 
(Alaghband-Zadeh et al. in preparation). We have zoomed in on the central 
~ 5000 kms -1 of the spectrum for clarity, although for completeness we 
also show the entire 2 GHz wide spectrum as an inset. 
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Nv and a continuum break (see Cop pin et al ] |2009l for the optical 
spectrum). Note that this apparent redshifting of Lya compared to 
the systemic redshift of the source (as measured from the nebular 
emission) is routinely seen in samples of 2 g z g 3 star-forming 
galaxies with offs ets of < AvLyq >^ 450 km s _1 compared to the 
systemic redshift dSteidel et alj|2010l) . The CO redshift also corre- 
sponds to the redshift of the [Oil] emission at z — 4.751 ± 0.005 
from a VLT spectrum (Alaghband-Zadeh et al. in preparation). This 
CO detection has provided definitive proof tha t the previously iden- 
tified counterpart at z = 4.76 for this SMG bv lCoppin et al.i ( 2009) 
is correct. 

Averaging the remaining channels in the ~ 1.6 GHz (129- 
channel) line-free region of the spectrum reveals no significant con- 
tinuum detection down to a 1 a sensitivity of 0.013 mJy. 

4 INFERRED GAS, DYNAMICAL, AND TOTAL 
BARYONIC MASS OF THE SYSTEM 

We calculate the line luminosity and the total cold gas mass 
(H2+He) from the int egrated CO(2-l) line flux following 
Solo mon & Vanden Boutl J2005t) in order to determine the amount 
of available fuel for the current star formation episode. We find a 
line luminosity of Lco( 2 -i) = (2.0 ± 0.4) x 10 10 Kkm s _1 pc 2 . 
We then assume both that the gas is thermalised for the J = 2- 

1 transition (i.e. assuming a constant line temperature brightness 
ratio, Lco(2-l)=Lco(l - 0)) and a CO-to-H2 conversion factor of 
a = O.8M0(Kkms -1 pc 2 )" 1 (which is appropriate for ultralumi- 
nous infrared galax ies and SMGs; e.g. iDownes & Solomon|[l998l ; 
iTacconi et al . 2008), yielding a total cold gas mass of M gas =( 1 .6 ± 
0.3) x 10 10 M o . This CO(2-l) gas mass is comparable to that 
derived for SMGs from higher-J transitions (typically J = 3— 

2 or J = 4-3). If LESS J033229.4 is a typical SMG, these CO 
J = 2-1 observations suggest that we are unlikely to be miss- 
ing m assive reservoirs of molecular gas aro und SMGs (see also 
e.g. ICarilli et alj|20ld: llyison et al.ll201ol ; cf. lHainline et afll2006l : 
Papadopou los et alj|2010h . 

We calculate a dynamical mass of the system based on the 
observ ed CO(2-l) line width following Solomon & Vande nBoutl 
d2005h : Mdyn sin 2 i=233.5A np WHM xR (M©), where i is the in- 



Figure 2. CO(2-l) intensity contours of the line peak in a 90 kms ^wide 
channel ATCA spectrum overlaid on a true-colour ground-based HAWK- 
I JK image of LESS J033229.4. The CO(2-l) solid line contours begin at 
2.5 a and increase in steps of 0.5 cr, where a = 0.14 mJy (the channel rms). 
Equivalent negative contours are shown by dotted curves. The CO source 
appears to be spatially unresolved with the ss 12 x 9 arcsec 2 beam (~ 80 X 
60 kpc FWHM; shown in the bottom right-hand corner). The CO centroid is 
consistent with the radio-identified optical V-band dropout counterpart of 
the SMG (encircled in green and featured as an inset in the top right-hand 
corner) within the % 2.3 arcsec positional uncertainty. 



clination of the gas disk, A i>fwhm is the CO FWHM (in km s" 1 ), 
and R is the radius of the CO emitting region (in pc). Since we do 
not have a direct constraint on the spatial extent of the CO emission 
from our unresolved detection, we assume that the CO traces the 
stellar light (as seems an appropriate assumption for z ~ 2 SMGs; 
ISwinbank etafl l201Ch . which has a seeing-corrected FWHM of 

0. 5" ± 0.2" in the HAWK-I AT-band imaging (or equivalently a 
FWHM ~ 4 kpc at z = 4.76). The best-fitting Gaussian FWHM 
of (160 ± 40) kms -1 implies a strict lower limit on the dynam- 
ical mass of M dyn (< 2 kpc) sin 2 i=( 1.2 ± 0.6) x 1O 1O M , as- 
suming a 2 kpc radius. Since we do not have any constraints on 

1, in this letter we assume i = 30° (appropriate for randomly 
inclined disks in a sample of galaxies) and explicitly include the 
i and R dependency in our derived value: Md yn (< 2 kpc) = 
(4.8 ± 2.4) x 10 10 (0.25/sin 2 i)(i?/2kpc)M Q . We caution that 
Mdyn carries significant uncertainties which will require higher- 
resolution observations to accurately determine these quantities. 
For instance, we note that the CO line FWHM is about a factor of 3 
narro wer than typically observed for z ~ 2 SMGs ( Gre ve et"al] 
2005), although there are examples of SMGs and high-redshift 
Lyman-break galaxies (LB Gs) with CO emi s sion as narrow as 
our observed FWHM (e.g. iFraver e7allll999l: iBaker et al .1120041; 
Coppin et al. 2007; Stanway et al. 2008). QSOs typically have nar- 
rower linewidths than SMGs by a factor of ~ 2-3, consistent with 
the optically-selected QSOs being more closely inclined to the sky 
plane (showing its central AGN more clearly) than typical SMGs 
(iGreve et alJioolJCoppin et al.l2008l ; lcarilli & Wand2006l) . Thus 
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the CO line narrowness of LESS J033229.4 (whose Hubble Space 
Telescope (HST) rest-frame UV compact morp hology and stron: 
Nv emission indicates that it hosts an AGN; ICoppin et al.11200 1 
compared with other SMGs could indicate that i ~ 10° . Assum- 
ing that LESS J033229.4 is more face-on, adopting i ~ 10°, yields 
M dyn (< 2kpc) = (4.0 ± 2.0) x IO^Mq. However, the line 
width differences could also indicate that LESS J033229.4 is at 
a later stage of a merger or has a different galaxy mass or size 
than typical z ~ 2 SMGs, which could be ex plored with high- 
resolution imaging of the gas distribution (e.g. iGreve et"aL 2005; 
Carilli & Wang 2006). If instead we assume that the gas config- 
uration is spherical (with a uniform distibution) rather than disk- 
like, then the implied enclosed dynamical mass would be Md yn (< 
2kpc) = (1.1 ±0.6) x 1O 1O M . 

The combination of the gas and stellar mass estimates for 
LESS J033229.4 can be used to calculate the total b aryonic mass 
of the system. The stellar mass o f the system fr o m ICoppin et al.l 
(2009) has b een reconfirmed by [ Wardlow et al.l ( 120100 to be ^ 
1 x 10 11 M m . IWardlow et alj ( 12010) has included the new HAWK- 
I J and -band photometry, and adjusted the derived mass for 
a ISalpeterl dl955t> IMF for compatibility with the IMF assumed 
in deriving the SFR in § 15.11 We caution that this quantity has 
an uncertainty of a factor of < 5 even before considering a po- 
tential AGN contribution to the rest-frame near-infrared emission 
and so should be considered an approximate upper limit. Taken to- 
gether, these estimates imply a total baryonic mass within 2kpc 



of M ba 



=M eas +M s 



1. 



1 x lO^M, 



which is consistent with the dynamical mass estimated above from 
the CO emission given the considerable uncertainties. 

5 DISCUSSION AND CONCLUSIONS 

5.1 z ~ 4-5 SMGs are higher redshift analogues of z ~ 2 
SMGs 

How representative is LESS J033229.4 of the few known z ~ 4- 
5 submillimetre-selected sources and the more abundant z ~ 
2 SMG population? Based on the photometric constraints for 
LESS J033229.4 which trac e an SEP consistent with local star for- 
mation dominated ULIRGs, Cop pin et alj |2009) derive a dust mass 
estimate of Ma ~ 5 x 10 8 Mq and a far-infrared luminosity of 
Lir ~ 6 x 10 12 L Q , implying a SFF0 of ~ 1000 M Q yr _1 , close 
to th e median luminosity of the z ~ 2 and z > 4 SMG populations 
(e.g. lKovacset1ul2006f) . When combined with our gas mass esti- 
mate above it follows that LESS J033229.4 has a star formation ef- 



ficiency (SFE) of « 250 L Q (K km s" 1 pc 2 



. a gas-to-dynamical 



mass fraction of /=Mg as /M, 



dyn ~ 0.3(sin i/0.25), and a gas-to- 
dust mass ratio of ~ 30. Assuming that the star formation follows 
the 7f-band light, we estimate a star formation surface density of 
Ssfr ~ 100 Mq yr _1 kpc~ 2 , which is similar to the intense cen- 
tral s tarburst mode inferred for SMGs at z ~ 2 (e.g. Sm ail et all 
120031) . 

The source properties of LESS J033229.4 derived from the 
CO line emission are very similar to those of the numerous SMGs 
at z ~ 2 «M gas >~ 3.0 x 10 10 M Q , <M dyn >~ 1.2 x 10 11 M , 
/ ~ 0.25, with SFEs of ~ 450 ± 170L B (Kkms" 1 pc 2 )" 1 
and M gas /M dust ~ 60: IGreve et al] 120051 ; iSovacs et alj l2006t) 
and also the rarer emerging high-redshift tail z > 4 SMGs 



3 The implied SFR is calcu lated following iKennicuf] jl998l) . assuming a 
100 Myr burst lifetime and a Salpeter 1 1955) initial mass function. 



dSchinnerer et al.l2008l ; lDaddi et alj|2009al lbl). Combining the sam- 
ples of CO-detected z > 4 SMGs, it follows that they lie within the 
scatter of the Lir-L' co relation from lGreve et al.1 fe005h for z ~ 2 
SMGs (assuming a constant line brightness ratio), with z ~ 4-5 
SMGs spanning a range of Lir ~ 0.6-3 x 10 13 L and L^o ~ 2- 
6 x 10 10 Kkms -1 pc 2 . It thus appears that within the z ~ 1-5 
SMG population, Lir/Lco is constant with redshift, indicating that 
SMGs at z ~ 4-5 are consistent with being higher-redshift ana- 
logues of SMGs at z ~ 2, forming stars with similar efficiencies. 

Overall this comparison suggests that SMGs at z ~ 2 and 
2 > 4 are equally evolved, have similar reservoirs of gas, similar 
star formation efficiencies, and similar fractions of baryons in cold 
gas as stars. 

5.2 The potential descendents of 4 < z < 5 SMGs 

LESS J033229.4 appears to be a luminous massive starburst, form- 
ing stars at a rate of ~ 1000 Mq yr~\ with a total baryonic mass 
of ~ 1 x 10 Mq, with properties representative of SMGs at 
2 ~ 2-4. The combination of its compact morphology and high 
SFR is potentially consistent with the small sizes of g 1-2 kpc 
claimed for a populatio n of massive extremely dense old quies- 
cent z ~ 2 galax i es (e.g.lDaddi et alj|2005l: iToft et alJ|2007l 120091 : 
IZirm et al.l [20071 : iBuitrago et alj 120081 : ICimatti et al.N2008l) . Such 
compact dense galaxies at z ~ 2 have implied formation redshifts 
of 2 ~ 4-5 and are consistent with being the descendents of a 
population of gas-rich, highly dissipative mergers, such as z > 4 
SMGs. 

A critical question we can now address is whether 
LESS J033229.4 (assumed to be representative of z > 4 SMGs) 
has the baryonic mass and gas consumption timescale necessary 
to be a 'prototypical' progenitor of the luminous 'red and dead' 
galaxy population found at z ~ 3? The baryonic content of 
LESS J033229.4 is roughly equivalent to the typical stellar mas s 
of a giant elliptical galaxy (1 x 10 11 M ; Marc hesini et al.ll2009T) . 
given the large uncertainties in some of our mass estimates. As- 
suming that the molecular gas reservoir is fueling the star forma- 
tion within LESS J033229.4, then it will have enough gas to sus- 
tain the current star formation episode for Tdepletion ~M gas /SFR~ 
1.6 x 10 10 M Q /1000M Q yr" 1 ~ 16Myr, assuming 100% ef- 
ficiency. We also compare repletion with the time to form the 
current stellar mass of the system. At the current SFR, we esti- 
mate a r formation ~M stars /SFR~ 1 x lO^MQ/lOOOMoyi'- 1 ~ 
100 Myr, which is comparable to the assu med age of th e stella r 
population and burst in the model used in ICoppin et al.1 (2009). 
Although the gas consumption timescale appears to be relatively 
short, it is sim ilar to the gas dep letion timescales of ~ 100 Myr 
for SMGs (e.g. lTacconi et al.ll2008l) . It follows that we are catching 
this SMG approximately half way through its current star formation 
episode, with the majority of the galaxy mass already in the form 
of stars, representative of a major stage in galaxy formation. It thus 
seems plausible that LESS J033229.4 would display the observed 
properties of being 'red and dead' by z ~ 3 since the current stel- 
lar population would age another > 1 Gyr between z — 4.76 and 
2 = 3, assuming that no further star formation occurs. 

5.3 Observed mass constraints of 4 < z < 5 SMGs on galaxy 
formation models 

The presence of even the small observed (and likely still incom- 
plete) volume density of the five z ~ 4-5 SMGs, 2 x 10~ 7 Mpc" 3 
(conservatively assuming a selection function that is uniformly 
sens itive across z = 4-8 due to the negative /("-correction; 
e.g. lBlainet"aUl2002h . poses a challenge to current hierarchical 
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models (see Cop pin et alj2009l) . However, the prediction of submm 
fluxes from the models depends on the assumed IMF and dust prop- 
erties in the galaxies, making a straightforward comparison be- 
tween the properties of the observed and model SMGs fundamen- 
tally problematic. A much more straightforward and reliable test is 
thus provided by comparisons of the total baryonic (col d gas+stars) 
or dy namical masses for these young galaxies (e.g. iGenzel et al.l 
2003), which we now investigate using observational Mbary and 
Mdyn constraints of the five CO-detected z > 4 SMGs, with a 
range of M bary ~ 1-3 x 10 11 M Q and M dy „ ~ 0.5-9 x 10 11 M . 
In Fig.[3]we compare the observational and theoretical cosmic vol- 
ume density of galaxies with Mb ary & lx 10 11 Mq (the lower end 
of the distribution among the sample), using high-redshift SMGs 
as examples of massive baryonic galaxies with a uniform selection 
function out to z ~ 8. This exe rcise demonstrates in a less-biased 
way that the lBaughetai] (2005 ) semianalytic GALFORM model is 
currently consistent with our current observed limits on the volume 
densities of M bary > 1 x 10 11 M SMGs at z > 4. 

These observations provide a view of a massive gas reservoir 
in a galaxy forming just 1 Gyr after the Big Bang, giving a sneak 
preview of the capabilities of ALMA for studying the formation 
and evolution of the earliest galaxies. However, statistically signifi- 
cant samples of spectroscopically-confirmed SMGs are required to 
determine the true incidence of SMGs at z > 4, and far-infrared 
and (sub)millimetre continuum surveys are first needed to identify 
large samples of SMGs that can become the targets of interfero- 
metric follow-up efforts such as those described here. Deep mul- 
tiwavelength near-infrared imaging and sensitive (preferably) low 
J-transition CO mapping campaigns are necessary for deriving ac- 
curate stellar and gas masses, respectively: the key to determining 
the contribution of SMGs to the baryonic mass of the Universe at 
z > 4. 
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